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FILMCAL 

AN  ESR  SIMULATION  PROGRAM  FOR  TRANSITION  METALS  IN 

PHTHALOCYANINE  FILMS 


M.  Dale  Pace  and  Stephanie  M.  Masumura 
Code  6120 

Naval  Research  Laboratory 
Washington,  D.  C.  20375-5000 


I .  INTRODUCTION 


The  applications  of  ESR  simulations  of  transition  metal 
ions  in  phthalocyanines  and  porphyrins  include  oxygen  transport 
in  hemoglobin,  and  semiconducting  films  for  microelectronic 
devices.  One  application  is  to  understand  chemiresistor  gas  sensors 
which  are  quartz  subtrates  5mm  by  8mm  with  4  0  interdigital  gold 
fingers  coated  with  Langmuir-Blodgett  films  of  a  semiconducting 
phthalocyanine  (See  figure  1.) 
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Figure  1 

Schematic  diagram  of  a  chemiresistor  with  forU' 
interdigital  gold  electrodes  coated  with  a  Langmuir- 
Blodgett  film. 

Manuscnpi  approved  November  21.  1988. 


These  devices  have  intrinsic  conductivities  which  change  when 
gaseous  vapors  are  absorbed  into  the  surface.  Each  gas  gives  a 
distinct  change  of  conductivity  for  a  given  film.  With  an  array 
of  such  detectors  placed  aboard  a  ship,  one  could  monitor  the 
levels  of  toxic  gases  present. 

Unfortunately,  the  reaction  time  of  the  films  is  from  5  to  25 
minutes  with  forty-five  layers  of  phthalocyanine.  In  order  for  a 
film  to  recover  its  initial  conductivity  hours  are  needed. 

One  factor  which  has  a  bearing  on  the  time  of  recovery  is  how  well 
ordered  the  phthalocyanine  coating  is  on  the  device.  The  method 
used  to  transfer  a  layer  of  phthalocyanine  onto  the  device  is  the 
Langmuir-Blodgett  technique.  This  method  transfers  a  layer  from 
a  water  surface  onto  the  substrate  (See  figure  2.) 
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Figure  2 

A  transference  of  the  phthalocyanine  coating  from  the 
water  surface  to  a  substrate.  The  !s  represent  polarized 


molecules. 

This  method  has  its  drawbacks  since  getting  the  coating  to  be 
perfectly  ordered  is  difficult.  This  is  due  to  changes  in  the 
surface  tension  of  the  water  which  must  be  monitored  at  all  times 
to  insure  that  the  polarized  coating  is  not  disturbed.  This 
computer  program  is  helpful,  for  it  allows  the  user  to  check  how 
well  ordered  the  film  coating  is  once  it  is  on  the  chemiresistor . 
This  is  due  to  the  fact  that  an  ordered  film's  spectrum  will  be 
completely  different  from  that  of  an  unordered  film.  The  program 
allows  the  user  to  check  possible  molecular  order  paramerers  by 
means  of  spectral  comparison  between  the  simulated  spectra  and  the 
experimental  spectra. 


Figure  3 

A  metal-subsdtuted  phthalocyanine,  where  the  M 
represents  a  metaJ  atom  and  N  a  nitrogen. 

Another  application  is  with  porphyrins  which  serve  as 
structures  of  hemoglobin,  chlorophyll,  and  some  enzymes. 
Hemoglobins  within  the  bloodstream  complex  with  oxygen  bound  along 
the  Z-axis  and  the  porphyrin  in  the  XY  plane.  When  oxygen  is 
unbound  the  iron  is  in  the  Fe^'^  state  having  five  d  electrons  with 
the  electrons  unpaired.  This  gives  a  paramagnetic  molecule  which 
is  detected  using  ESR.  Accurate  simulation  of  the  spectrum  and 
comparison  to  the  experimental  spectrum  can  be  used  to  determine 
the  quantity  of  oxygen  within  the  measured  blood  sample.  So  it  is 
clear  that  this  program  has  many  practical  uses. 


II.  THEORY 


The  program  was  written  to  calculate  spectra  for  semiordered 
systems  having  S=l/2,  or  in  the  case  of  S=1  calculating  the  lines 
but  not  the  intensities  (The  S=1  case  does  not  formally  calculate 
the  lineshapes  but  it  does  calculate  the  magnetic  fields  at  which 
the  transitions  occur  for  a  completely  random  spectrum.  Future 
versions  of  the  program  will  include  the  lineshape) .  For  the  S=l/2 
case  the  following  spin  Hamiltonian  is  used. 

H=g(3H'S  -H  S'A^'I^  +  S‘A^'1^ 


where  S=l/2;  M  denotes  different  metal  ions  having 
different  hyperfine  values  (A)  and  different  nuclear  spin 
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numbers  1=1/2,  3/2,  5/2,  or  7/2  and  denotes  the 

superhyper fine  coupling  to  the  nitrogen  ligands  having 
1^=1.  N  has  values  of  1  to  4  denoting  the  four  near-neighbor 
nitrogen  ligands.  The  nuclear  Zeeman  interaction  and  the 

guadrupole  interaction  are  not  included.  Thescs  values  are 
relatively  small  and  in  most  cases  do  not  affect  the  spectral 
lineshapes. 


The  program  begins  by  setting  up  the  parameters  for  the 
spectral  calculations.  The  first  question  prompts  the  user  to 
choose  between  several  options  available,  the  first  being  whether 
or  not  to  include  the  spin=l  case.  This  S=1  case  is  used  for 
situations  where  dimerization  of  the  molecules  results  in  spin 
pairing  so  that  each  paramagnetic  species  has  S=l.  The  spectrum 
produced  is  derived  by  solving  the  following 
matrix  for  the  value  of  D  and  E. 


(D/3)  +  gPHn 
g{3H(l4-tm) 

E 


g3H(l-im) 

-2D/3 

.ggHg+im) 

i2 


E 

(D/3)  -  gPHn 


To  formally  calculate  the  transition  probabilities  a  polynomial 
equation  is  solved  for  the  field  values  at  which  pairs  of  energy 
levels  are  separated  by  g  H.  These  fields  are  inserted  into  the 
above  matrix  which  must  be  diagonalized  to  obtain  the  state 
energies  and  vectors.  The  transition  intensities  are  computed  from 
these  equations.  This  must  be  carried  out  for  each  value  of 
magnetic  field  orientation.  This  is  time  consuming  and  for  this 
program  the  explicit  calculation  of  the  transition  intensities  has 
been  omitted.  Instead,  we  have  programmed  the  general  solution  of 
the  field  values  of  the  transitions  when  D  and  E  are  input.  This 
can  be  used  to  aid  in  assignment  of  a  completely  random  S=1  system 
when  the  experimental  spectrum  is  available. 

The  case  most  often  used  is  S=l/2.  Within  a  normal 
phthalocyanine  the  metal  atom  in  the  center  is  the  site  of  unpaired 
electron  density.  For  S=l/2  this  results  in  a  hyperfine 

interaction  between  the  nucleus  of  the  metal  and  the  unpaired 
electron  spin.  Depending  on  the  bonding  and  electronegativities 
of  the  nitrogens  the  unpaired  spin  density  can  be  delocalized 
sufficiently  so  that  there  is  an  interaction  with  the  four 
surrounding  nitrogen  nuclei.  These  four  nitrogen  interactions 
add  eight  extra  lines  to  the  spectrum  and  must  be  accounted  for  in 
the  lineshape  procedure. 


The  choice  of  including  second-order  perturbation  is  the 
next  option.  When  there  is  a  large  hyperfine  interaction  as  with 
the  transition  metal  ions  the  lines  of  the  spectrum  will  be  shifted 
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slightly  from  their  first-order  positions.  To  first  order  the 
distances  between  the  transitions  are  equal.  However,  in  the  case 
of  second  order  the  transitions  are  shifted  causing  the  lines  to 
become  unevenly  spaced.  This  effect  is  shown  in  figure  4.  without 
the  second  order,  metals  like  vanadyl,  which  has  a  nuclear  spin 
moment  of  7/2  have  computed  spectra  which  are  incorrect. 

After  making  these  selections,  the  user  is  requested  to  enter 
the  number  of  nuclear  spin  isotopes.  Then,  for  each  isotope,  the 
nuclear  moment  and  %  abundance  is  needed.  When  there  is  only  one 
isotope,  any  constant  is  acceptable  for  the  natural  abundance. 

The  sweep  width  which  is  the  range  of  points  displayed  on  the 
X-axis  is  then  needed.  This  entry  is  in  gauss.  The  program  will 
function  with  any  sweep  width  size  although  it  begins  to  lose 
accuracy  if  it  is  not  in  the  range  of  100  to  10000.  The  linewidth, 
the  next  piece  of  data,  is  the  distance  from  peak-to-peak  (in 
gauss)  of  a  line. 

The  next  data  is  the  g^^,  g  ,  and  g  grouping.  These  values 
are  used  for  the  calculation^of  the  location  of  the  hyperfine  lines 
and  are  generally  near  the  value  of  2.  The  lower  this  value,  the 
larger  the  gauss  index;  in  other  words,  the  spectrum  is  shifted 
farther  to  the  right.  The  A^,  and  K  values  are  next  input  in 
gauss.  These  are  the  principal  values  or  the  hyperfine  interaction 
tensor. 

Then,  the  option  is  selected  and  the  superhyperf  ine 

couplings  are  input.  These  are  referred  to  as  A^^,  A  ,  and  A^ 
and  have  the  same  meaning  as  described  above.  ^ 

Once  these  data  have  been  entered  the  program  asks  for  angle 
rotations.  This  permits  calculations  of  different  orientations 
of  the  molecule  within  3-D  space.  This  space  is  referenced  by  the 
X,  Y,  and  Z  axes,  with  the  angles  of  tilt  from  each  being  theta 
(0) /  phi  ($) ,  and  rho  (P) ,  respectively. 


These  choices  are  important  because  a  spectrum  of  an  ordered 
sample  will  look  different  depending  on  the  angle  at  which  you  scan 
when  using  the  ESR  instrument.  If  a  powder  spectrum  is  being 
computed,  the  program  should  be  allowed  to  scan  all  possible  angle 
combinations.  The  reason  is  that  the  polymer  sample  is  disordered 
and  all  orientations  of  the  polymer  molecules  must  be  accounted 
for.  The  spectra  are  then  averaged  together  to  achieve  the  final 
result. 

For  Langmuir-Blodgett  films,  the  coating  is  very  layered 
and  ordered,  and  this  must  be  accounted  for  when  computing  the 
spectrum.  The  only  problem  with  this  is  that  the  molecules  usually 
align  themselves  not  completely  perpendicular  or  parallel 
to  the  surface,  but  rather  at  a  bit  of  an  angle.  Sigma  is  the 
standard  deviation  of  the  most  probable  tilt  angle.  The  option  of 
entering  the  most  probable  angle  allows  the  user  to  account  for 
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FIGURE  ^ 
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this.  If  the  angle  is  unknown  test  spectra  can  be  computed  and 
compared  to  the  experimental  spectra  recorded  at  different 
orientations  of  the  film  with  respect  to  the  applied  magnetic  field 
direction.  By  matching  the  experimental  and  computed  spectra  at 
all  orientations,  an  assignment  of  the  most  probable  tilt  angle 
and  standard  deviation  can  be  made.  This  assumes  that  one  has 
accurately  kept  account  of  the  coordinate  axes  of  the  experimental 
reference  frame. 

After  programming  these  options,  the  program  requires  input 
for  graphing.  We  decided  to  use  the  DISSPLA  graphics  package  that 
was  on  the  system.  This  graphing  package  uses  calls  to  a  calcomp 
plotter  to  do  the  plotting. 


III.  EXAMPLES  OF  SIMULATIONS 

FILMCAL  is  an  ESP  (electron  spin  resonance)  spectral 
simulation  program  for  calculating  semiordered  systems  having 
S=l/2  and  1=1/2,  3/2,  5/2,  or  7/2.  The  program  is  specifically 
written  to  calculate  transition  metal  ion  ESR  spectra  of 
phthalocyanines  and  porphyrins  which  have  a  central  metal  ion 
coordinated  to  four  near-neighbor  atoms  of  the  molecular  rings. 
The  hallmark  of  this  program  is  the  flexibility  to  calculate 
different  orientational  distributions  of  the  molecular  ri.ngs 
(phthalocyanines  and  porphyrins  are  planar  ring-like  molecules) 
with  respect  to  a  fixed  XYZ  frame  of  reference  representing  a 
substrate  surface.  The  program  has  been  extensively  used  to 
simulate  ESR  spectra  of  thin  organometallic  films  composed  of 
metal-substituted  phthalocyanines;  thereby  determining  the 
molecular  packing  arrangement  of  the  molecules  which  comprise  these 
films.  This  gives  a  direct  and  simple  procedure  for  monitoring  the 
molecular  order  of  thin  films  by  using  the  ESR  spectra. 

The  program  requires  minimal  input.  An  XYZ  frame  of  reference 
is  assigned  as  the  substrate  coordinate  system  with  X  normal  to  the 
substrate  surface  which  is  taken  as  the  XY  plane.  A  second 
coordinate  system  abc  is  assigned  as  the  molecular  axis  system  of 
the  phthalocyanine  molecule  with  c  normal  to  the  molecular  plane 
which  is  the  ab  plane.  This  is  shown  schematically  in  figure  5. 
An  angle  rho  (P) ,  defined  as  the  tilt  angle,  is  the  angle  between 
the  Z  axis  of  the  substrate  frame  and  the  c  axis  of  the  molecular 
frame  assuming  that  the  X  and  a  axes  always  lie  coincident  along 
the  same  direction.  Only  rotation  about  the  X  and  a  axes  direction 
is  allowed.  This  provides  a  means  of  specifying  different  tilt 
angles  of  the  molecular  plane  (ab)  wit):  respect  to  the  substrate 
plane  (XY) ,  Direction  cosines  are  used  to  specify  the  direction 
of  the  polarizing  magnetic  field  (Hq)  with  respect  to  the  XYZ 
reference  frame.  Polar  coordinates  are  used  where  theta  (0)  is  the 
angle  between  a  component  vector  of  Hq  in  the  ZX  plane  and  the  Z 
axis.  Phi  ($)  is  the  angle  of  a  component  vector  of  Hq  in  the  XY 
plane  from  the  X  axis.  This  permits  Hq  to  be  specified  in  any 
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direction  and  is  analogous  to  an  experiment  where  the  substrate  is 
rotated  with  respect  to  Hq. 

This  gives  sufficient  degrees  of  freedom  to  specify  any 
orientation  of  the  molecular  reference  system  with  respect  to 
XYZ  and  Hq.  For  example,  to  calculate  a  spectrum  corresponding 
to  the  molecular  system  coincident  with  the  substrate  system  and 
with  Hq  normal  to  the  substrate  plane  (ie.  along  Z)  we  input: 


A.  EXAMPLE  OF  A  POLYCRYSTALLINE  SPECTRUM  OF  VO^^ 

(SEE  SPECTRUM  FIG.  6) 

Is  this  an  S^l  case?  0  [where  l=yes  and  0=no] 

Include  N-14  interaction?  0 
Include  2nd  order?  1 
Number  of  isotopes?  1 

Spin  1,  nuclear  moment=1.0  [when  only  one  isotope  any  constant 

is  acceptable] 

Spin  number  =3.5  [VO  nuclear  spin  of  7/2] 

%  abundance®  100. 

Sweep  width®  2000. 

Linewidth®  10.0 
=  1.9951 
Gy  =  1.9951 
G^  =  1.9673 
Metal  ion  HFC  values 

=  59.7  [input  is  in  gauss] 

Ay  =  59.7 
A^  =  170.1 

-  Theta  - 

Start  =  0 
Stop  =  90 
Step  =  2 
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^entaaquo  Vanadyl(IV)  ion  in  FVozen  Aqueous  Solution 
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2639.97  2609.94  3139.91  3380. 8&  3639.84  3889.81  4139.78  4389.75 

MAGNETIC  FIELD  (GAUSS) 

Figure  6. 


Calibration  =  l  [Calibration  allows  for  steps  between  integer 

angles.  For  example  an  input  might  be: 

Start=o 
Stop  =  1000 
Step  =  2 

Calibration  =  0.09 

This  would  step  from  0  to  90  degrees 
by  500  steps  ie.  1000  x  0.09  =  90.] 

- Phi - 


Pstart  =  0 
Pstop  =90 
Pstep  =  2 
Calibration  =1 

-  Angle  of  Tilt  From  2  Axis 


Start  =  0 
Stop  =90 
Step  =  2 
Calibration  =  l 

-  Most  Probable  Angle  of  Tilt  - 

Tilt  angle  =  0  [If  Sigma  is  larger  than  100  any  value  of  the  most 

probable  tilt  angle  is  okay.  This  will  weight 
all  tilt  angles  equally] 

Sigma  =  1000 
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2390.00  2639.97  2889.94  3139.91  3389.88  3639.84  3889.8 


B.  EXAMPLE  OF  A  POLY CRYSTALLINE  SPECTRUM  OF  Cu^^ 

(SEE  SPECTRUM  FIG.  7) 

This  example  computes  a  polycrystalline  ESR  spectrum  of 
copper  phthalocyanine. 

Is  this  a  S=1  case?  0 

Include  N-14?  1 

Include  2nd  order?  0 

Number  of  isotopes?  2 

Spin  1,  nuclear  moment  =  1.588 

Spin  nvimber  =  1,5 

%  abundance  =  69.2 

Sweep  width  =  2000. 

Linewidth  =  20. 

Gx  =  2.045 
Gy  =  2.045 
G^  =  2.187 

Metal  ion  HFC  values 
=  35.0 
Ay  =  35.0 
A^j  =  228.3 
N-14  SHFC  values 
=  17.1 
Ay  =  17.1 
A^  =  15.5 

-  Theta  - 

Start  =  0 
Stop  =  90 
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step  =  2 

Calibration  =  1 


- Phi - 

Start  =  0 
Stop  =  90 
Step  =  2 
Calibration  *  1 

-  Angle  of  Tilt  from  Z  Axis  - 

Start  =  0 
Stop  =  90 

Step  =  2 
Calibration  =  1 

-  Most  Probable  Angle  of  Tilt 

Tilt  angle  =  0 
Sigma  =  1000. 


C.  EXAMPLE  OF  LANGMUIR-BLODGETT  FILM  SPECTRA 
(SEE  SPECTRA  FIG.  8) 

Disordered  ESR  spectra  of  an  LB  film  of  VOPCXi,  were  calculated 
for  two  alignments  of  Hq  (parallel  and  perpendicular  to  the 
surface  of  the  film) .  For  these  calculations  the  most  probable 
tilt  angle  (p)  is  held  constant  at  80*.  The  tilt  angle  refers  to 
the  angle  between  the  plane  of  the  phthalocyanine  ring  and  the 
planar  surface  of  the  substrate.  The  standard  deviation 
is  varied  to  simulate  the  disorder.  The  numbers  10’,  20’,  30’, 
50’  refer  to  the  number  of  degrees  of  deviation  of  tilt  angle  from 
the  value  p.  The  disorder  increases  reading  down  each  column. 
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I V .  Listing  of  ESR  spectral  simulation  program 


PROGRAM  FILMCAL 

C - - - 

C  Filmcal  -  An  ESR  Simulation  Program 
C  Developed  by  Dale  Pace  &  Stephanie  Masumura 
C  Naval  Research  Laboratory,  July  1988 
C 

C  Filmcal  is  an  ESR  simulation  program  which  calculates 
C  spectra  of  transition  metal  ions  having  S=l/2  and  1=  1/2, 

C  3/2,  5/2,  or  7/2.  The  program  was  developed  in  order  to 
C  simulate  ESR  spectra  of  organometallic  films  in  which 
C  transition  metal  ions  act  as  spin  probes.  The  program 
C  operates  by  summing  spectra  of  different  values  of  theta 
C  (the  angle  between  the  polarizing  magnetic  field  and  the 
C  z  principal  axis  of  the  g-  and  hfc  tensors  in  the  zx  plane) 

C  and  phi  (the  angle  between  Ho  and  the  x  principal  axis  in 
C  the  xy  plane)  and  rho  (the  tilt  angle  of  the  g-  and  hfc 

C  tensors  from  the  xyz  reference  frame).  A  tilt  angle  of  the 
C  g-value  and  hfc  tensors  from  the  xyz  principal  axis  refer- 
C  ence  frame  is  specified  by  the  angle  rho  in  order  to  simulate 
C  different  orientations  of  molecules  which  make  up  the  film. 

C  A  weighting  factor,  W,  is  ^cified  as  a  gaussian  function 
C  in  order  to  give  different  dikributions  of  the  molecular 
C  orientations  within  the  film. 

C  An  approximate  simulation  of  5=1  spectra  is  included  as 
C  an  option.  This  simulation  calculates  field  positions  of 
C  spectra  having  different  D  and  E  values  to  assist  with  S=1 
C  assignments.  The  formal  lineshapes  and  relative  transition 
C  intensities  are  not  included  in  this  version  of  the  program. 

C 

C - 

REAL  LINE,LINE1, FLAG ,TERM1,TERM2MP. WINDOW 
REAL  G,HFC,THETA,H,SC ALE,M7(8000). V( 3 ,3 ,3 ).T( 3 .3 .3 ) 

REAL  PLANCK,FREQ,PSI,RHO,CAL,PCAL,TCAL,W.SIGMA 
REAL  SRHO,CRHO,SSRHO,SCRHO,SPSI,CPSI.SSPSI,SCPSl 
REAL  AB(3),M0(3).SPIN(3)4NWF(9)MRATI0.M(8) 

REAL  XARAY(8()00),HIGHNUM,LOWNUM.LW.HR.XINCRE 
REAL  DEN.NUMJL\NGE.INCREJv18(8000).HZ,NUMBER.HX{8) 
INTEGER  LD^E2.START,ST0P.STEP.PSTEP.PSTART.PST0P.TSTEP 
INTEGER  RHOl  JlH02,n,JJ.LT0TAL,N,SPAN 
C 

CON^/=9.3377E-5 

PLANCK=6.6252E-27 

FREQ=9.5E9 

BETA=0.92731E-20 

GE=2.0023 

C - INITIALIZATIONS - 

DO  101=1,3 
DO  11  J=l,3 
DO  12K=1.3 
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V(U4C)=0 

12  CONTINUE 
11  CONTINUE 
10  CONTINUE 

00  441=1,8000 
M7(I)=0 
M8a)=0 
XARAYa)=I 
44  CONTINUE 
PSI=0.0 

: - SET-UP - 

PRINT  *, 'SET-UP' 

PRINT  *,TES=1  NO=0' 

PRINT  '(A,$)','  IS  THIS  A  S=1  CASE?  ' 

ACCEPT  *.FS 
IF  (FS.EQ.l)THEN 
PRINT  '(A,$)'.'  D=  ' 

ACCEPT  *,D 
PRINT  '(A,$)','  E=  ' 

ACCEPT  *,E 

PRINT  '(a’$)','  SWEEPWIDTH=  ' 

ACCEPT  *,SW 

PRINT  LINEWIDTH=  ' 

ACCEPT  *,LW 
WINDOW=3390-(SW/2) 

GOTO  9999 
ENDBF 

PRINT  '(A,$)','  INCLUDE  N-14  INTERACTION?  ' 

ACCEPT  *,F1 

PRINT  '(A,$)','  INCLUDE  2ND  ORDER?  ' 

ACCEPT  *,F2 

: - ISOTOPE/GRAPHING  DATA - 

PRINT  '(A,$)';  NUMBER  OF  NUCLEAR  SPIN  ISOTOPES? 
ACCEPT  *,N 
DO  13  1=1. N 
PRINT  2,1 

2  FORMATC  SPIN’ ,13,'  NUCLEAR  MOMENT=  ',$) 

ACCEPT  *MO(I) 

PRINT  '(A,$)',’  SPIN  NUMBER=(DECIMAL)  ’ 

ACCEPT  *.SPIN(I) 

PRINT  '(A,$)','  %  ABUNDANCE=  ' 

ACCEPT  *  AB(I) 

13  CONTINUE 

PRINT  '(A.$)','  SWEEP  WIDTH=  ' 

ACCEPT  * SW 
SCALE=NIOT(8000/SW) 

WINDOW=3390-{SW/2) 

PRINT’(A.$v;  LrNEWIDTH=  ’ 

ACCEPT  *.LW 

: - DATA  VALUES - 

9999  PRINT  '(A,$)'.'  GX=  ' 

ACCEPT  *.V(1, 1,1) 

PRINT  ’(A,$)','  GY=  ' 

ACCEPT  *.V(2.2.1) 


PRINT  GZ=  ' 

ACCEPT  *.V(3,3,1) 

IF  (FS.EQ.l)  GOTO  9998 

C - HFC  DATA - 

PRINT  *, 'METAL  ION  HFC  VALUES' 
PRINT  '(A,$)','  AX=  ’ 

ACCEPT  ♦,V(  1,1, 2) 

PRINT  ’(A,$)’,’  AY=  ’ 

ACCEPT  *,V(2,2,2) 

PRINT  ’(A,$)','  AZ=  ■ 

ACCEPT  *,V(3,3,2) 

C - N-14  DATA - 

IF  (Fl.EQ.l)  THEN 
PRINT  *,N-14  SHFC  VALUES' 
PRINT  '(A,$)','  AX=  ■ 

ACCEPT  *,V(  1,1, 3) 

PRINT  ’(A,$)',’  AY=  ’ 

ACCEPT  *,V(2,2,3) 

PRINT  '(A,$)','  AZ=  ' 

ACCEPT  *,V(3,3,3) 

NWF(1)=1 

NWF(2)=4 

NWF(3)=10 

NWF(4)=16 

NWF(5)=19 

NWF(6)=16 

NWF(7>=10 

NWF(8)=4 

NWF(9)=1 

ENDIF 


DO  25  K=l,3 
DO  26  1=1,3 
DO  27  J=1.3 

V(LJ,K)  =  V(I,J,K)  *  V(I,J,K) 

PRINT  UJ JC,V(1,J,K) 

1  FORMAT  {’  V(M2,',M2,',’.I2,')=’.G16.8) 
27  CONTINUE 
26  CONTINUE 
25  CONTINUE 

C - ANGLES - 

PRINT  - THETA - ' 

PRINT  '(A,$)'.'  START=  ' 

ACCEPT  ♦.START 
PRINT  '(A,$)'.'  STOP=  ■ 

ACCEPT  ♦.STOP 
PRINT  '(A,$)’,'  STEP=  ’ 

ACrFPT  *  STEP 

PRINT  '(A.'S)'.'  CALIBRA'i  10N=  ' 

.-ACCEPT  ♦, CAL 

PRINT  ♦.' - PHI - ' 

PRINT  '(A,$)'.'  PSTART=  ' 

ACCEPT  ♦.  PSTART 
PRINT  '(A,$)'.'  PSTOP=  ' 
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ACCEPT  *^STOP 
PRINT  '(A,$)',’  PSTEP=  ' 

ArCFPT  *  PSTFP 

PRINT  •(A,$)';  CALIBRATION=  ’ 

ACCEPT  *,PCAL 

PRINT  - ’ 

PRINT  *  .’ANGLE  OF  TILT  FROM  Z- AXIS' 

PRINT  '(A,$)',’  START=  ’ 

ACCEPT  *IIH01 
PRINT  ’(A,$)’,'  STOP=  ’ 

ACCEPT  *.RH02 
PRINT  ’(A,$)’.’  STEP=  ' 

ArCFPT  *  TSTFP 

PRINT  '(A,’$)’,'  CALIBRATION=  ’ 

ACCEPT  *.TCAL 

PRINT  - ’ 

PRINT  ♦,'MOST  PROBABLE  TILT 
PRINT  ’(A,$)','  TILT  ANGLE=  ’ 

ACCEPT  *MP 
MP  =  MP*  .01745 
PRINT  '(A,$)’,’  SIGMA=  ' 

ACCEPT  *,SIGMA 
SIGMA=SIGMA*SIGMA 
DO  1000  I=START, STOP  .STEP 
THETA  =  .01745  *  I  *  CAL 
DO  2000  J=PSTARTJ>STOP4’STEP 
pm  =  .01745  *  J  *  PCAL 
X=SIN(THETA)*COS(Pm) 

Y=SIN(THETA  )*SIN(Pm) 

Z=COS(THETA) 

DO  50  n=RHOLRH02.TSTEP 
RHO  =  -.01745  *  n  *  TCAL 
W=(MP+RHO)/.01745 
W=EXP(-(W  *  W)  /  (2  *  SIGMA)) 

C - CALCULATIONS  WITHOUT  THE  SECOND  ORDER - 

IF  (F2.EQ.0)  THEN 
SRHO=SIN(RHO) 

CRHO=COS(RHO) 

SSRHO=SIN(RHO)*SIN(RHO) 

SCRHO=COS(RHO)*COS(RHO) 

SPSI=SIN(PSI) 

CPSI=COS(PSI) 

SSPSI=SIN(PSI)*SIN(PSI) 

SCPSI=COS(PSI)*COS(PSI) 

DO  70  K=l,3 

C - INTERMEDIATE  CALCULATIONS  FOR  G  AND  HFC 

T(  1 .1  JC)=(  V(  1 , 1  JK:)*SCPSI)+((SSPSI)*  •  ( V(2,2  jC)*SSRHO 
1  +(V(3.3JC)*SCRH0)))) 

T(l,2jC)=SPSI*((V(2,2X:)*SRHO*CRHO)-(V(3.3.K)*SRHO*CRHO)) 
T(1.3,K)=(-V(1,1  jC)*CPSI*SPSI)+<CPSI*SPSP((V(2,2JC)*SSRHO) 

1  +{V(3.34C)*SCRHO))) 

C 

T(2,l.K)=(SPSI*((V(2.2.K)*CRHO*SRHO-(V(3,3.K) 

1  *SRHO*CRHO)))) 
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T(2,24C)=(  V(2,2  JC)*  SCRHO)-K  V(3 ,3  .K)*SSRHO ) 
T(2,3JC)=CPSI*((V(2,2,K)*CRH0*SRH0HV(3,3.K)*SRH0*CRH0)) 
C 

T(3 ,1 4?^)=((- V(  1 , 1  ,K)*SPSI*CPSI)-KSPSI*CPSI 
1  *((V(2,2,K)*SSRH0)+{V(3,3,K)*SCRH0)))) 

T(3.24C)=CPSI*((V(2,2,KrSRHO*C:RHO)-(V(3,3.K)*SRHO*CRHO)) 

T(3,3JC)=(V(l,l,K)*SSPSI)+{CPSI*((V(2,2,K)*SSRHO) 

1  -KV(3,3JC)*SCRHO))) 

70  CONTINUE 

C - CALCULATION  OF  G  AND  HFC  WTTH  VALUES  FROM  ABOVE - 

K=1 

G=SQRT(X*X*T(1,1^)  +  X*Y*T(2,14C)  +  X*Z*T(3,14C) 

1  +  X*  Y*T(  1 ,2,K)  +  Y*  Y*T(2,2^)  +  Y*Z*T(3 ,2,K) 

1  +  X*Z*T(1,3  JC)  +  Y*Z*T(2,3,K)  +  Z*Z*T(3.3^)) 

C 

K=2 

HFC=SQRT(X*X*T(1,1  JC)  +  X*Y*T(2,1  JC)  +  X*Z*T(3,1.K) 

1  +  X*Y*T(1,2^)  +  Y*Y*T(2,24C)  +  Z*Y*T(3,2,K) 

1  +  X*Z*T(1,3^)  +  Y*Z*T(2,3.K)  +  Z*Z*T(3,3  JC)) 

C - N-14  INTERACTION - 

IF(FLEQ.1)THEN 

K=3 

NHFC=SQRT(X*X*T(1,1,K)+X*Y*T(2,LK)+X*Z*T(3,1,K) 

I  +X*  Y*T(  1  ,2,K)+Y*  Y*T(2,2  JtC)+Z*  Y*T(3 ,2,K) 

I  +Z*X*T(1,3.K)+Z*Y*T(2,34^>+Z*Z*T(3,3JC)) 

ENDIF 

H=PLANCK*FREQ/(G*BETA) 

D0  3000in=l.N 
LTOTAL=2*SPIN(ni)+ 1 
LT0TAL=LT0TAL-1 
MRATIO=MO(ni)/MO(  I ) 

LLOW=(HFC*LTOTAL*MRATIO)/2 
DO  210  JJJ=0,LTOTAL 
TERMl=HFC*JJJ*MRATIO 
LINE=H-LLOW+TERM  1 

C - N-14  INTERACTION  PLACEMENT  OF  NITROGEN  LINES - 

IF(F1.EQ.1)THEN 
DO  211  KKK=l,9 
TERM2=NHFC*(KKK-5) 

LINE  1 =LINE-t-TERM2 
LINE2=NINT((LrNEl-WINDOW)*SCALE) 
IF(LINE2.LE.8000.AND.LINE2.GT.0)  M7(LrNE2)=M7(LINE2) 

1  -)-W*AB(m)*NWF(KKK) 

211  CONTINUE 

C - PLACEMENT  OF  UNES/NON-NTTROGEN - 

ELSE 

LINE2=NINT((LINE-WINDOW)*SCALE) 

IF  (LINE2.LE.8000)  M7(LrNE2)=M7(LrNE2)+W*.\B(in) 

ENDIF 

C - 

210  CONTINUE 
3000  CONTINUE 
C - SECOND  ORDER  CASE - 
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on 


ELSE 

CTHETA=SIN(THETA)*SlN(PHI)*SIN(RHO)+COS(THETA)*COS(RHO) 
STHETA=ABS(  l-(CTHETA*CrHETA)) 

STHETA=SQRT(STHFrA) 

G=SQRT(  V(3 ,3 . 1  )*CTHETA*CTHET  A+ V(  1.1,1  )*  STHETA*  STHETA ) 
HFC=(SQRT((V(3 ,3 ,2)*  V(3 ,3 , 1  )*CTHETA*CTHET  A) 

1  +(V(1,1,2)*V(1,1,1)*STHETA*STHETA))) 

HFC=HFC/G 

C - SECOND  ORDERyN-14  INTERACTION - 

IF(F1.EQ.1)THEN 

NHFC=SQRT((V(3 .3 .3)*  V(3,3, 1)*CTHETA*CTHETA) 

1  +(V(1,1,3)*V(1,1,1)*STHETA*STHETA)) 

ENDEF 

C - 

DO  60  ni-l.N 
LTOTAL=2*SPIN(ni)+ 1 
LTOTAL=LTOTAL-l 
MRATIO=MO(in)/MO(  1 ) 

MI=-(SPIN(in)) 

DO  61  JJJ=0,LTOTAL 

A=(  V(  1 ,1 ,2)*CONV*CONV*(V(3,3^)*CONV*CONV 
1  +HFC*HFC*CONV*CONV)) 

A=A/(4*3. 169*HFC*HFC*CONV*CONV) 

A=A*(15.75.(MI*MI)) 

B=((V(l,l,2)*CONV*CONV-V(3,3,2)*CONV*CONV)/(HFC*CONV)) 
B=B*B 

B=(.5*B)/.3169 

E==STHETA*STHETA*CTHETA*CrHETA 
B=B*E*MI*MI 

LINE=(0.3i69+HFC*MI*CONV*MRATIO-A-B)/(G*4.6686E-5) 

- SECOND  ORDER/N-14  INTERACTION - 

- PLACEMENT  OF  NITROGEN  LINES - 

IF(F1.EQ.1)THEN 
DO  62  KKK=1.9 
TERM2=NHFC*(KKK-5) 

LINE  1 =LINE+TERM2 
LINE2=NINT((LINE  1  -  WINDOW)*SCALE) 

IF(LINE2.LE.8000)  M7(LINE2)=M7(LINE2)+W*AB(m)*N^TlKKK) 

62  CONTINUE 

C - PLACEMENT  OF  LINES  W/0  N-14 - 

ELSE 

LINE2=NINT(ABS((LINE-WINDOW)*SCALE)) 

IF  (LINE2.LE.8000.AND.LINE2.GT.O)  M7(LINE2)=M7(LINE2)+W*.\B(in) 
ENDIF 

C - 

\n=Nn+i 

61  CONTINUE 
60  CONTINUE 
ENDIF 

C - END  OF  SECOND  ORDER - 

50  CONTINUE 
2000  CONTINUE 
1000  CONTINL^ 

GO  TO  3010 
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c - s=l  CASE - 

9998  BETA2=4.6686E-5 

HO=(PLANCK*FREQ)/{GE*BErA) 

DP=D/(GE*BETA2) 

EP=E/(GE*BCTA2) 

IF  ((DP**2-(3*(EP**2))).GT.O)  DPS=SQRT(DP*-^2-(3*(EP**2))) 
SCALE  =NINT(8000/SW) 
HX(l)=(GE/V(3,3,l))**2*(((HO-DP)**2HEP*EP)) 
HX(2)=(GE/V(2,2,1))**2*((H0-DP-EP)*(H0-(EP*2))) 

HX(3)=(GE/V(  1,1,1  ))**2*((HO-DP+EP)*(HO+{EP*2) ) ) 

HX(4)=(GE/V(  1,1,1  ))**2*((HO+DP-EP)*(HO-(EP*2) )) 
HX(5)=(GE/V(2,2,1))**2*((H0+DP+EP)*(H0+(EP*2))) 
HX(6)=(GE/V(3,3,l))**2*((HO+DP)**2-(EP*EP)) 
HX(7)=((HO/2)**2-(DPS**2/3)) 

HX(8)=(HO**2-(DPS**2/3)) 

DO  JJJ=1,8 

IF  (HX(JJJ).GT.O.AND.HX(JJJ).NE.O)  THEN 
HX(JJJ)=SQRT(HX(JJJ)) 

SINDEX=NINT(((HX(JJJ)-WINDOW)/((8000/SCALE)/8000))- 1 ) 
IF(SINDEX.GT.O.AND.SINDEX.LT.(WINDOW+SW)) 

1  M7(SINDEX)=M7(SINDEX)+1 
ENDIF 
ENDDO 

C - LINE  SHAPE  DETERMINATION - 

3010  PRINT  *,'DETERMINING  THE  LINE  SHAPE  OF  THE  SPECTRUM' 
SPAN=2*SCALE*LW 
LW=(LW/SW)*8000 
DO  4000  1=1,8000 
IF(M7(I).GT.O)  THEN 
DO  5000  J=0,SPAN 
HZ=XARAYa+J)*SCALE 
HR=XARAY(I)*SCALE 
NUM=2*LW*  LW*(  HZ-HR ) 

DEN=LW*LW+((HZ-HR)*(HZ-HR)) 

DEN=DEN*DEN 

IF  (J+I.LE.8000)  M8(I+J)=-M7a)*(NUM/DEN)+M8(I+J) 

IF  (J.GT.l)THEN 

IF  (I-J.GT.O)  M8(I-J)=M7(I)*(NUM/DEN)+M8(I-J) 

ENDIF 

5000  CONTINUE 
ENDIF 

4000  CONTINUE 

C - DETERMINATION  OF  EXTREME  Y  V.ALLES - 

HIGHNUM=M8(1) 

LOWNUM=M8(l) 

DO  33  JJJ=2,8000 
IF  (M8(JJJ).GT.HIGHNLTvI)  THEN 
HIGHNUM=M8(JJJ) 

ELSE  IF  (M8(JJJ).LT.LOWNUM)  THEN 
LOWNUM=M8(JJJ) 

ENDIF 

33  CONTINUE 

C - SCALING  LT  OR  DOWN  OF  VALUES - 

C - TO  CREATE  A  BETTER  GRAPH - 


NUMB  ER= ABS(LO  WNUM) 

IF  (fflGHNUM.GT.NUMBER)  THEN 
NUMBER=HIGHNUM 
ENDIF 

IF  (NUMBER.LT.60.AND.NUMBER.NE.0)  THEN 
999  DO  36  KKK=  1.8000 
M8(KKK)=M8(KKK)*  1 0 
36  CONTINUE 

NUMBER=NUMBER*  10 
IF  (NUMBER.LT.60)  GOTO  999 
ELSE  IF  (NUMBER.GT.1000)  THEN 
199  DO  39  KKK=  1,8000 
M8(KKK)=M8(KKK)/1 0 
39  CONTINUE 

NUMB  ER=NUMB  ER/1 0 
IF  (NUMBER.GT.1000)  GOTO  199 
ENDIF 

C - DETERMINATION  OF  A  SUITABLE  RANGE  FOR  GRAPHING - 

RANGE=(INT(NUMB  ER/1  OOH 1  )*  1 00*  1 . 

C - DETERMINATION  OF  A  SUITABLE  INCREMENT  FOR  Y-AXIS - 

INCRE=100. 

NUM=RANGE/100 
IF  (NUM  MOD  2.EQ.0)THEN 
INCRE=200. 

ENDIF 

C - SCALES  THE  X  ARRAY  CORRECTLY  INTO  GAUSS - 

RATE=((8000/SCALE)/8000) 

XARAY(l)=WINDOW 
DO  49  m=2,8000 

XARAY(in)=XARAY(in- 1  )-i-RATE 
49  CONTINUE 

XINCRE=(XARAY(8000)-WINDOWV8 

C - GRAPHING  OF  THE  SPECTRA  WITH  DISSPLA  9,2 - 

PRINT  'PREPARING  TO  GRAPH  THE  SPECTRUM' 

CALL  COMPRS 
CALL  PAGE(11.,8.5) 

CALLTRIPLX 

CALL  TITLE  (TITLESMOO.'MAGNETIC  HELD  (GAUSS)$',100, 

1  'SIGNAL  INTENSITY$'.100.8..5.5) 

CALL  GRAF  (WINDOW ,XINCREJCARAY(8000),-(RANGE),INCR£JLAMGE) 
CALL  CLTIVE(XARAYM8.8000,0) 

CALL  ENDPL(O) 

CALL  DONEPL 
END 
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VI.  DISKETTE 


The  source  code  FILMCAL.FOR  compiles  using  MICROSOFT  FORTRAN 
using  the  FL  /Gt  command.  The  field  values  of  array  XARAY(8000) 
and  the  intensities  of  array  M8(8000)  are  output  in  SYS$OUT. 
These  arrays  may  be  dimensioned  to  different  lengths 
by  changing  the  approriate  dimension  statements  and  loop 
counters  in  the  program.  A  plotting  program  such  as  ASYST  must 
be  used  to  plot  the  SYS$OUT  file  and  is  not  included  in  the 
program  code. 

An  diskette  of  the  program  is  available  from: 

M.  D.  Pace 
Code  6120 

Naval  Research  Laboratory 
4555  Overlook  Ave. ,  S.  W. 

Washington,  D.  C.  20375-5000 
(202) -767-3239 
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